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Abstract. 

We present high spectral resolution optical spectra of three young intermediate mass stars, in all of which we spectrally resolve 
the 6300 A [OI] emission line. Two of these have a double peaked line profile. We fit these data with a simple model of the 
[OI] emission caused by photo-dissociation of OH molecules in the upper layer of a circumstellar disk by stellar UV radiation 
and thus translate the Doppler broadened [OI] emission profile into an amount of emission as a function of distance from 
the central star. The resulting spectra are in agreement with the expected disk shapes as derived from their spectral energy 
distribution. We find evidence for shadowing by an inner rim in the disk surrounding HD101412 and see a flaring disk structure 
in HD179218 while the [OI] spectrum of HD135344 B is more complex. The [OI] emission starts for all three targets at 
velocities corresponding to their dust sublimation radius and extends up to radii of 10-90 AU. This shows that this method 
can be a valuable tool in the future investigation of circumstellar disks. 



1. Introduction 

HD101412, HD135344 B and HD179218 are three isolated 
young intermediate mass stars. Two of them, HD101412 and 
HD 1792 18 are Herbig Ae/Be (HAEBE) stars while the third, 
HD135344 B (often incorrectly referred to as HD135344 in the 
literature), is a classical T Tauri star of spectral type F. HAEBE 
stars are the more massive counterparts of T Tauri pre-main- 
sequence stars with typical masses of 2 - 8 M . Their spectral 
energy distribution (SED) is characterized by the presence of 
an infrared excess due to thermal re-emission of circumstel- 
lar dust which is thought to be the signature of a circumstel- 
lar disk. Previous studies of these stars and their circumstellar 



99% of the disk mass, the gas. To check this coupling and 
to get a better picture of circumstellar disks in general, it is 
necessary to study the dust as well as the gas located in cir- 
cumstellar disks simultaniously. Next to the abundance, there 
are other reasons such as its potential to probe the dynam- 
ics and the physical and chemical structure of disks that illus- 
trate why the circumstellar gas is of interest. Historically it has 
been very difficult to study the circumstellar gas, but the de- 
velopment of bigger mirrors and better spectrographs in the 
recent years has opened the way for seve ral gas diagnostics . 
These include hot and cold CO gas (e.g. iBik & Thil d2004: 
iBlum et al.1 (12004 : iBlake & Boogertl d2004 and iNaiita et al 
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disks have focused mainly on this dust (e.g Hillenbrand et al 
dl992h. iNatta et al l d200ll) , iMannings & Sargent! d 19971) and 
lMeeusetalJd200ll) ). 

Dust however contains only 1 % of the mass of the disk 
and it is unclear whether this dust is coupled to the other 



(12003b ). hot water (e.g.lCarret alJ d2004) and m o lecula r hy 
drogen (e.g. ICarmona et al.l d2008l) : iBitneretal. (2007) and 



Martin-ZaidietajJ <|2007)). 



Acke et all d2005l) and lAcke & van den Ancked d2006l) 
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have added molecular oxygen to this list. They have developed 
a method for examining excited neutral oxygen atoms in the 
surface of a flared, rotating passive disk that they have tested on 
2 HAEBE stars: HD97048 and HD100546. With their method, 
they resolve the disk rotation and the distribution of the emit- 
ting gas in the surface layer of the circumstellar disk and also 
find evidence for a gap in the disk around HD 100546 that they 
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interpret as due to the presence of an orbiting body of planetary Table 1. Log of spectroscopic observations with UVES 
mass within the circumstellar disk. 



In this paper we apply the method of lAcke et al.l d2005l) on 
three more stars, one of which has a non-flaring (flat) disk ge- 
ometry. We first describe these stars, together with our obser- 
vations and data reduction in Sections [2] and [3] then proceed 
to describe the method and the resulting spectra in Section [4] 
and discuss these results in Section [5] We have also observed 
our targets with VLT/MIDI to get spatially resolved observa- 
tions of the 10 pm silicate bump whose emission originates 
from warm dust in the same region as the [OI] emission. The 
results of these interferometric observations will be presented 
in Fedele et al. (2008, in prep; Paper II). 



2. Observations and Data Reduction 

We obtained high spectral resolution echelle spectra of 
HD 1 1 4 1 2 , HD 1 35 344 B and HD 1 792 1 8 during 4 nights with 
the KUE YEN 8.2 meter ESQ VLT coupled to the UVES spec- 
trograph (Dek ker et al.l l2000). The observation log is presented 
in table Q] The data was reduced with the UVEsQ pipelined 
which performs wavelength calibration, order extraction, back- 
ground subtraction and flatfield correction. The observations 
consist of 4 data sets per target. Half of these sets were taken 
with the slit unrotated (0 degrees), and the other half with 
the slit rotated 90 degrees. The spectra have a spectral reso- 
lution of — 77000, determined using the telluric absorp- 
tion lines around the 6300 A [OI] feature. All three sources 
are spatially unresolved. We derive upper limits on their an- 
gular size using a 3 cr detecti o n limi t and followi n g the same 



Hirthetal. (1997J) and 



procedure as e.g. iHirfhetalJ dl994h . 
Acke & van den Anckerl d2006l) . The upper limits to the size 
of the emitting region are 0.25 arcsec corresponding to 39 AU 
for HD101412, 0.26 arcsec or 36 AU for HD135344 B and 
0.22 arcsec or 53 AU for HD179218 . These sizes represent 
the radii that contain 95% of the flux in the 6300 A [OI] line. 
The observations were made at two different dates to coincide 
with VTLI/MIDI observations of the same targets in the mid- 
infrared (N-band). Since the spectra taken with different slit 
orientation show negligible variation they have been averaged. 
The data procured at different times however do show temporal 
change and therefore are handled independently. They will be 
referred to as Bl (May 15 - 17) and B2 (June 16 & July 3, 4) 
in this paper. 

We determined the radial velocity for our targets from the 
difference between the measured line center and the theoretical 
line center for 9 to 3 1 atmospheric absorption lines observed in 
the same spectrum. The errors have been determined by boot- 
strapping the data set, i.e. by estimating the variability of the 
radial velocity from the variability of this velocity between sub 
samples. The size of the sub samples used are 90% of the to- 
tal sample and we have estimated the variability by calculating 
the mean of the subsample a 1000 times. The results are shown 



Object 


Date 


UT 


Exp. time 


PA 


B 


S/N 




dd/mm/yy 


hh:mm:ss 


[s] 


[°] 






HD101412 


15/05/06 


23:37:48 


300 





Bl 


174 




15/05/06 


23:43:45 


300 





Bl 


199 




15/05/06 


23:50:04 


300 





Bl 


189 




15/05/06 


23:56:01 


300 





Bl 


175 




17/05/06 


00:39:17 


300 


90 


Bl 


183 




17/05/06 


00:45:15 


300 


90 


Bl 


207 




17/05/06 


00:51:13 


300 


90 


Bl 


155 




17/05/06 


00:57:09 


300 


90 


Bl 


199 




03/07/06 


23:52:01 


300 





B2 


149 




03/07/06 


23:58:00 


300 





B2 


141 




04/07/06 


00:03:58 


300 





B2 


161 




04/07/06 


00:09:53 


300 





B2 


151 




04/07/06 


00:18:58 


300 


90 


B2 


137 




04/07/06 


00:24:53 


300 


90 


B2 


144 




04/07/06 


00:30:48 


300 


90 


B2 


135 




04/07/06 


00:36:46 


300 


90 


B2 


155 


HD 135344 B 


17/05/06 


04:19:43 


300 





Bl 


197 




17/05/06 


04:25:39 


300 





Bl 


196 




17/05/06 


04:35:08 


300 


90 


Bl 


234 




17/05/06 


04:41:04 


300 


90 


Bl 


262 




16/06/06 


03:49:40 


300 





B2 


144 




16/06/06 


03:55:36 


300 





B2 


143 




16/06/06 


04:06:12 


300 


90 


B2 


227 




16/06/06 


04:12:29 


300 


90 


B2 


240 


HD179218 


17/05/06 


06:50:00 


180 





Bl 


309 




17/05/06 


06:56:28 


180 


90 


Bl 


237 




16/06/06 


04:22:43 


180 





B2 


293 




16/06/06 


04:28:05 


180 


90 


B2 


251 



1 http://www.eso.org/instruments/uves/ 

2 http://www.eso.org/projects/dfs/dfs-shared/web/uves/uves-pipe- 
recipes.html 



in Table [2] The [OI] 6300 A line is located in between sev- 
eral telluric absorption lines. To remove these, we fit and fill 
the absorption lines with a Gaussian, and average the telluric 
line region over 3/2 FWHM to suppress telluric line residuals. 
These telluric absorption lines only overlap the [OI] spectral 
region on the red wing of HD 10 14 12. The results are shown in 
Figures [T]-|3] In these Figures we also show the location of the 
6300 A [OI] sky emission line. This line has been removed by 
calculating the sky emission from averaging above and below 
the echelle order and subtracting this from the final spectrum. 
Because the sky emission is very weak compared to the flux 
collected from the targets, the procedure has no significant ef- 
fect on either the resulting spectrum or the signal to noise. 

The spectra were also corrected for their radial and 
barycentric velocity. Then we normalize the spectra and con- 
vert them to velocity profiles around the central wavelength 
of the [OI]6300 A line (6300.304 A). The resulting spectra of 
HD101412 and HD179218 are shown in Figures |4]and|6]while 
the spectrum of HD 135344 B in Figure|5]needs to undergo one 
more reduction step to correct for the underlying photospheric 
absorption lines as discussed in section l4.2.2l 

The spectra of HD101412 and HD179218 both show a dou- 
ble peaked [OI] line profile indicating the presence of a circum- 
stellar disk and temporal variance indicating disk inhomogene- 
ity. HD 10 14 12 shows a maximum signal of 10 % above the 
continuum, has 60 km s _1 wide wings and displays next to the 



G. van der Plas et al.: The Structure of Protoplanetary Disks Surrounding Three Young Intermediate Mass Stars 



3 



8000 : 



7000 : 



6000 r 



5000 



4000 



3000 E 





6296 6298 6300 6302 

wavelength[A] 



6304 




6298 6300 6302 

wavelength[A] 



6304 



Fig. 1. The spectrum for HD101412. The clipped region is 
shown in black for the Bl (upper line) and the B2 (lower 
line) spectra, and is also shown with horizontal bars above the 
clipped parts. The grey area denotes the location of the (sub- 
tracted) 6300 A [OI] sky line. The emission is not symmetric 
around this line because these spectra are not yet corrected for 
the radial and barycentric velocity 
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Fig. 2. Same as FigureQ]for HD135344 B . Note that the 6300 
A [OI] line is mixed with broad atmospheric absorption lines. 
The over plotted horizontal dotted line indicates the continuum 
and the telluric absorption lines of the B 1 data have been trun- 
cated for so it wont overlap with the B2 data. 



Fig.3. Same as Figure Q]for HD179218 , the telluric absorp- 
tion lines of the Bl data have been truncated for so they wont 
overlap with the B2 data. 
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Fig. 4. In the top window the velocity spectrum of HD101412 
centered around the [OI] line, after correcting for barycen- 
tric and radial velocity and normalization, for Bl (black) and 
B2(grey). The location of the clipped out telluric lines is shown 
by horizontal bars. The difference between both observations is 
shown in the bottom window. The interstellar absorption fea- 
ture at 50 km s _1 discussed in section|2]is clearly visible and re- 
placed with the blue shifted part of the same spectrum as shown 
with the dotted line. 



double peak two extra 'shoulders'. The data also show tempo- 
ral change with a more pronounced double peak in the B 1 data 
as well as more flux in the shoulders. We deduce that the dif- 
ference between the two lines on the red wing at 50 km s _1 in 
the B2 data seen in Figure |4] is not due to temporal variabil- 
ity because it drops below the continuum. It is also not present 
in the [OI] emission at 6363 A and thus real absorption. This 
line coincides with the telluric absorption in the Bl data and 
therefore does not show in that spectum. We have compared 
this spectrum with that of HD34364, another B9.5V star but 
failed to identify the absorption feauture, thus ruling out photo- 
spheric absorption. We also exclude telluric emission because 



the width of the line of 15 km s is larger than the telluric line 
width of 4 km s _1 . We therefore conclude that the nature of this 
absorption line is probably interstellar, but we cannot identify 
the exact origin. For further analysis we replace the affected 
part of the spectrum between 35 and 65 km s _1 with the -65 to 
-35 km s _1 data of the same spectrurr0. We show this replace- 
ment with the grey dotted line, while the original absorption is 
shown with a grey solid line in Figure [4] 

3 In order to compare how much this changes our final results we 
show the results obtained with the two seperate wings in Figure QT] 
as well as the results obtained with the above mentioned changes in 
Figure [T3l 
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Fig. 5. Same as Figureg]for HD 135344 B 
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Fig. 6. Same as Figure g] for HD 1792 18 



The [OI] emission line seen in HD 1792 18 rises a maximum 
of 5 % above the continuum, is 50 km s _1 wide and also shows 
intrinsic temporal change. The red wing of the line changes be- 
tween both data sets; This is not caused by normalizing effects. 

The [OI] emission in HD 135344 B has a maximum of ap- 
proximately 8 % above the continuum and a single peaked pro- 
file. Because of the overlapping photospheric absorption lines 
it is difficult to accurately determine the velocity broadening of 
the line. We estimate it to be inbetween 60 and 100 km s _1 . The 
[OI] line shows no descernible temporal change and, contrary 
to HD101412 and HD179218, peaks at 3.9 km s 1 (3 <x) away 
from the expected stellar photospheric zero velocity. 

3. Stellar Parameters 

The Spectral Energy Distributions (SED) of HAEBE stars can 
be divided into two groups based on the s hape of the mid-IR 



(20-100 jim) emission; groups I and II (IMeeus et al.l 120011) 



Group I sources have a rising mid-IR flux excess indicative 
of a flared disk, while group II sources display a more mod- 
est IR excess indicating a geometrically flat disk. This division 
between groups is made more quantitative by a classification 



a -9 




HD 101412 
Sp. = AOIIIe 
E(B-V) = 0.20 



log A H 

Fig. 7. The SED of HD 101412 compiled using literature pho- 
tometry. The solid line is a reddened Kurucz stellar atmosphere 
model fitted to the photometry of the central star. Also shown 
are SPITZER-IRS data from Bouwman et al. (2007, in prep.) 



using the LnjrILj r vs. IRAS [12] - [60] colour diagram (c.f. 
Acke et al. (2004)). We show the SED for each of our sources 



in Figures|7]-|9] By using this classification and comparing the 
SEDs, HD179218 and HD135344 B are group I sources. In the 
LnirIL, r vs. IRAS [12] - [60] diagram HD101412 is located 
close to the line separating the two groups. Although it is a 
group II source, it may be considered a transitional object in 
betwe en both groups. 

A lKuraczl(ll99ll) model atmosphere corresponding to the 
photospheric parameters of the central star is fitted through the 
data to represent the photospherical contribution. From this, the 
shape and excess emission from the disk become clear. The 
shortest wavelength for which the excess is discernible is 0.9 
fim for HD101412, 1.5 //m for HD135344 B and 1.0 //m for 
HD179218 . From here on, the excess emission for HD101412 
shows a steady decline that 'flares up' at 100 fim. This behavior 
is displayed by more IRAS sources at long wavelengths and is 
probably caused by the inclusion of interstellar material heated 
up by, but not directly associated with, the central star in the 
rather large (1 arcmin) IRAS beam. The other two SEDs show 
a re-brightening at mid-IR wavelengths. The fractional lumi- 
nosity of the dust compared to the stellar luminosity, Lire/X*, 
is listed in Table [2] 

To determine the projected rotational velocity, the v sin i, 
we have fitted a Doppler profile to mimic the broadening of the 
lines by stellar rotation to the rotational line profile of several 
atmospheric absorption lines. 



3.1. HD101412 

HD101412 has 



between B9.5Ve 



spectral type 

dHouk & Cowlevl Il975l) and AO III/IVe dGuimaraes et al 



2006) and has a distance of 160 pc. By integrating under 
the SED shown in Figure [7] (using a reddened Kurucz model 
to interpolate between observed data points in the optical, 
a spline in the infrared and a grey body to extrapolate they 
last data points), and adopting this distance of 160 pc, we 
derive a bolometric luminosity of 25 L Q for HD101412. 
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Name 


Spectral Type 


logT eff 


log L bol 


Lire 


M 


Distance 


Vrad 


V sin i 






log [K] 


log [L Q ] 


[£*] 


[M Q ] 


[pc] 


[kms" 1 ] 


[km s- 1 ] 


HD101412 


AOIIIe 


4.02 


1.40 


0.27 


2.3 ± 0.2 


160^ 


16.9± 0.2 


8 ± 1 


HD 135 344 B 


F4Ve 


3.82 


1.01 


0.43 


1.7 ± 0.2 


14Q* 


1.6+ 1.3 


75 + 5 


HD179218 


B9e 


4.02 


1.88 


0.28 


2.7 ± 0.3 


240^ 


15.4 ±2.3 


72 ± 5 



Name 


[OI] EW 


L([OI]) 


[OI] inner radius 


[OI] outer radius 


Inclination 


Group 




[A] 


[L e ] 


[AU] 


[AU] 


n 




HD101412 


-0.13± 0.01 


(9.5 ±0.5) 10- 5 


0.15 


10 


30 


II 


HD 135344 B 


-0.13 ±0.06 


(1.3 ±0.5) 10~ 4 


0.1 


90 


45f 


1 


HD179218 


-0.053± 0.005 


(3.4 ±0.3) 10- 4 


0.4 


65 


401 


I 





AckeetalJj2005l) 


b 


van Boekel et alj J20051) 


c 


van den Ancker et al. 1 19981) 


d 


Doucet et all 120061) 


c 


DentetalJj2005b 





HD 179218 




Sp. = B9e 




E(B-V) = 0.17 








\ 







log \ [fj.m] 



Fig. 8. Same as Figu re[7]for HD 1792 18 but the ISO -SWS spec- 
trum is shown here dAcke & van den Ancketfl2004 



authors also note the presence of emission bands due to the 
presence of crystalline silicates at 19, 24, 28 and 34 yum. This 
presence of strong PAH emission is, just like the presence 
of strong [OI] emission, unex pected for a group II source 
( Acke & van den Anckej2004 ). This reinforces our suggestion 
that the disk of HD101412 is in transition between flaring and 
being self shadowed. 

We estimate the mass of HD101412, derived by compar- 
ing the position in the HR di agram to the theoretical pre-main- 
sequence tracks of IPalla & Stahled d 19931) . to be 2.3 + 0.2 M Q . 
Gui maraes et al. ( 2006h find a v sin i of 7 + 1 km s _I by mod- 



eling high resolution (R = 48000) spectra, suggesting we see 
HD101412 almost pole-on. This agrees well with the value de- 
rived from our data of 8 + 1 km s _1 . 




HD 135344 
Sp. = F4Ve 
E(B-V) = 0.10 



log X [/im] 

Fig. 9. Same as Figure|7]for HD 135344 B 



The Spitzer-IRS spectrum of HD 101412 has an unusual 
shape due to the presence of strong PAH emission bands on 
top of a plateau-like 10 fim silicate feature over a weakly 
rising continuum (Bouwman et al. 2007, in prep.). The same 



3.2. HD 135344 B 

HD135344 B is a wide visual binary with a separation of 19.6 
arcsec located at a distance of 140 pc. In this paper we will 
discuss component HD 135344 B which has spectral type F4Ve 
and an estimated mass of 1.7 + 0.2 M . It has a bolometric lu- 
minosity of 10 L Q , obtained by fitting the SED shown in Figure 
[9] Bouwman et al. (2007, in prep.) note only very weak PAH 
features in the Spitzer-IRS spectrum of HD135344 B. They 
also note that its SED is quite unusual, wit h an IR excess start- 
ing at 1.5 yum, but a deep bump at 13 ^m. iDoucet et al. I d2006h 
have imaged HD 135344 B in the mid-infrared (20.5 //m) and 
find a disk with an inclination of 45° and a PA of 100 +10°. 
Based on these spatially resolve d observations they estimate 
the disk to extend up to 200 AU. iBrown et alJ ( 2007h classify 
the disk around HD135344 B as a 'cold disk', interpreting the 
lack of mid-IR emission in the SED as missing warm dust. 
From modeling the SED they derive the disk to have an inner 
rim starting at 0.18 AU, a gap (i.e. missing warm dust) from 
0.45 to 45 AU and outside that a gas rich disk. 
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A simple model fit to t he CO (J=3-2) double peaked line 
profile done by iDent et al.l ([2005 ) indicates an inner radius of 
the emitting gas < 10 AU, an outer radius of 75 + 5 AU, an esti- 
mated circumstellar dust mass of 1O~ 4 M and an inclination of 
11+2°. The value of this i nclina tion is notably different from 
the result of Douc et et al. I d2006l) and implies that the star ro 



tates a t breakup velocity. We will use the result of iDoucet et al 



(2006) in our model because it is obtained in a more direct way. 



3.3. HD179218 

HD 1792 18 is a B9e star at a distance of 240 pc and is the 
Herb ig star with the 2 nd large st known percentage of crystalline 
dust (Ivan Boekel et al. I l2005h . It has a bolometric luminosity of 
76 Lq, obtained with the SED shown in Figure[8] HD179218 is 
the star with the greatest richness in terms of number of infrared 
spectral features in the overview of infrared spectra of Herbig 
Ae/Be stars by Bouwman et al. (2007, in prep.). Its Spitzer-IRS 
spectrum is dominated by strong emission features due to the 
presence of abundant crystalline silicates. Polycyclic Aromatic 
Hydrocarbons are also det ected, in agreement with its cla ssifi- 
cation as a group I source ( Acke & van den Anc ker 2004). We 
estimate the mass based on the position in the HR diagram to 
be 2.7 + 0.3 M^. lGuimaraes et alj d2006l) find a v sin / of 72 + 
3 km s _1 agreeing well with the value derived from our data of 
72 + Skms- 1 . 



Elia et al.l d2004 fit the continuum spectrum of HD 1792 18 



with a model including a polar cavity, suggesting a disk-like 
geometry in the innermost part of the envelope. Moreover, it 
is noteworthy that their model corresponds to a dust absorp- 
tion coefficient exponent /? = 0.6 whi ch is typical of larg e dust 
grains and agrees with the results of lAcke et al.l d2004l) . This 
is in agreement with the idea that the grain size is expected 
to increase with time in circumstellar envelopes, as a result 
of the ongoin g growth proce sses activated in high-density en- 



Liu et al. I d2007l) have resolved warm dust using 



vironments. 

10 micron nulling interferometric observations around the star 
with a diameter of 27 + 5 AU. Their observations suggest cir- 
cular symmetry and thus a small inclination, but inclinations up 
to 45° are still within the error b ars. A simple m odel fit to the 
CO(J=3-2) line profile done by IDent et al.l (12005b indicates an 
outer radius of 120 + 20 AU, an inclination of 40+10° and an 
estimated circumstellar dust mass of 10~ 4 M o . 



4. Data Analysis 

4.1. Method 

In pioneering work on t he mechanism r espon sible for [OI] 
emission in young stars, lHartigan et al.l (L1995) modeled the 
broad [OI] 6300 A line emission seen in T Tauri stars as be- 
ing due to a combination of de nse stellar jets and a disk wind 
or magnetic accretion columns. lAcke et alj d2005l) have devel- 
oped a model that explains the narrow [OI] emission lines seen 
around many HAEBE stars. This simplified model assumes that 
the entire UV luminosity (2-13.6 eV) of the central star (repre- 
sented by a Kurucz model for the stellar photosphere) is used 
to photo dissociate OH in the upper layers of the disk. The disk 
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Fig. 10. Total mass of gas needed to explain the received [OI] 
6300 A flux vs. the relative strength of the 5577 A vs. 6300 
A line (bottom) in case of thermal emission. The correspond- 
ing temperature that is needed to explain this ratio is plot- 
ted on the upper x-axis. Also shown are the data points for 
the 3 stars (squares) and the function connecting the ratio of 
Equivalent Widths to the total gas mass for HD101412 (solid 
line). Because both lines are only resolved for HD101412 we 
show lower limits for the other two targets. The disks around 
our three targets would need to be unrealistically massive to ex- 
plain the observed [OI] emission as being due to thermal emis- 
sion from the disk 

structure (i.e. densities and temperatures) is computed with a 
Chiang & Goldreichl d 19971) type disk, with small olivine dust 
grains as the opacity source. The OH to H abundance is fixed a 
priori and constant at each location in the disk. The absorption 
of UV photons by the dust in the disk is taken into account 
in the same way as in the models of Storzer & Ho llenbach 
(2000) for photo dissociation regions along the direction of 
the in falling photons. The line profiles are deduced assum- 
ing Keplerian rotation and integration over the spatial coordi- 
n ates. Details about the method we employ here can be found 
in 



Acke et al. (2005) and I Acke & van den Anck er (2006). 



To verify the validity of this assumption for the three sys- 
tems we study here we compare the Line Flux (LF) of the 5577 
A [OI] emission to that of the 6300 A [OI] emission. If this 
emission is caused by thermal emission originating from a disk 
with a temperature of 500 K, the Boltzmann factor renders the 
ratio of these t wo very low (Figu re [TU 



We follow lAcke et alj d2005l) and lAcke & van den Ancker 
(2006) in concluding that the observed [OI] emission is caused 
by illumination of the upper disk layer by UV radiation. 

We expect that the pattern and intensity of the [OI] emis- 
sion is influenced by differences in disk geometry. The part of 
the disk that lies in the shadow of the puffed up inner rim char- 
acteristic for group II sources will emit less compared to the 
parts that are directly illuminated by the central star. Therefore 
we predict that there is a drop in intensity at the location of 
the shadow cast by the inner rim in HD 10 14 12. How extended 
this shadow is depends on when -and if- the disk rises out of 
the shadow again. We also expect the flaring disks to appeal" 
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larger in [OI] emission, because they are illuminated by stellar 
UV flux at greater distance and with more intensity due to their 
shape than the self shadowing disks. 

We note that the radius scale of our results is dependent on 
the inclination of the disk used in the model. Taking a smaller 
inclination angle shifts the emission region closer to the star. 
If we compare the onset of the [OI] emission with the location 
of the dust sublimation radius, we see that these agree for the 
inclinations used for HD 1 35344 B and HD 1 792 1 8 . We use this 
to guess an inclination for HD101412. Values between 24° and 
41° allow the inner rim to be at the dust sublimation radius. 
In the rest of this paper we will use a value of 30°. This low 
inclination is in agreement with both the low v sin i and the 
double peaked line profile seen in [OI]. 

4.2. Results 

4.2.1. HD101412 and HD179218 

By comparing the blue- and red shifted Doppler profiles, it is 
possible to distinguish between the parts of the disk that rotate 
towards us or away and thus track possible inhomogeneities. 
We have done this for HD 10 14 12 and HD 1792 18 and show the 
results in FiguresQT]and[T2]plotted together with the dust sub- 
limation radius for temperatures between 1200 and 2000 K in 
grey. The blue- and red-shifted wings of HD101412 look ap- 
proximately the same in the upper (Bl) and bottom (B2) graph 
apart from some excess flux at larger distances in the blue 
shifted part of the disk. HD 1792 18 looks mildly asymmetric, 
showing for the red shifted part emission that starts at a closer 
distance to the star and lacking intensity at larger distances 
compared to the blue shifted part. The blue side also seems to 
harbor more flux and is detectable up to larger distances. The 
results discussed below for HD101412 and HD179218 are ob- 
tained by averaging the red- and blue shifted profiles of both 
stars. 

The distribution of emitting atomic oxygen gas as a func- 
tion of radius for HD 10 14 12 is plotted in Figure [T3l and for 
HD 1792 18 in Figure [14] We have plotted the log - log depen- 
dency of the received flux as well as a semi log scale plot where 
we have normalized the intensity with the surface of the ring at 
distance R : S(R) where S(R) = n* (R 2 out - /??) with R out and 
Ri„ the outer and inner radius of the ring at distance R. In this 
sense I{R)xS{R) is the total luminosity in a ring at distance 
R and the fluxes of the inner and outer part of the disk can 
be directly compared this way. The error bars drawn in these 
graphs are solely a function of the signal to noise of the data. 
The [OI] emission of HD101412 starts within the expected lo- 
cation of the inner rim at 0. 15 AU and extends to 10 AU with a 
clear drop of almost 50 % in the I(R) vs. S(R) curve around 0.5 
AU. Since HD101412 is a group II source, characterized by a 
puffed up inner rim, and because the [OI] emission mechanism 
is directly proportional to (stellar) UV flux, the drop in inten- 
sity can be caused by the shadowing effect of the inner rim. 
Once the disk emerges out of the shadow (around 1 AU), the 
disk is again directly illuminated by the central star and keeps 
receiving enough flux up to 10 AU to be detectable. If we com- 
pare the [OI] emission of HD101412 with that expected from a 
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Fig. 11. The [OI] emission of the blue shifted (grey, dia- 
monds) and red shifted (black, triangles) halves of the disk of 
HD101412 . The two panels show the data taken at different 
times. The grey area is the suspected location of the inner rim. 
The horizontal dashed lines mark in the upper plot the location 
of the removed telluric absorption line (Figure Q]), and in the 
lower plot the part of the B2 spectrum that has been replaced 
as described in section|2]and Figure|4] 
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Fig. 12. Same as FigureQUfor HD179218 

flared disk in Keplerian rotation, shown with light grey crosses 
in FigureQj] one sees that for a flared disk the emission should 
have been detected up to much larger radii. We interpret these 
discrepancies with the flared disk model that both can be ex- 
plained by a self-shadowed disk; the local drop in intensity at 
0.6 AU and the dissapearance of emission at 10 AU, as direct 
evidence that in the disk surrounding HD101412, the inner rim 
indeed shadows part of the outer disk - as expected for a group 
II source. 

The intensity vs. radius distribution for HD 1792 18, a group 
I source, looks different. It starts at 0.3 AU from the central 
star, extends^up to 65 AU and has no clearly discernible local 

4 This number given here is larger than the in section [2] inferred 
upper limit for the disk size. This is caused by the fact that the disk 
emission is concentrated towards the center and drops of exponentially 
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Fig. 13. The intensity versus radius graph of HD101412 show- 
ing the distribution of [OI] emission as a function of distance 
from the central star. The two lines represent data taken on dif- 
ferent dates; Bl is plotted in black with squares and B2 in grey 
with diamonds. The grey area is the suspected location of the 
inner rim and the light grey crossed line represents a simple 
model for [OI] emission from a flared disk in Keplerian rota- 
tion. In the upper plot we show the intensity normalized to the 
disk surface in a semi-log plot, as to better compare the inner 
and outer parts of the disk. A drop in intensity can clearly be 
seen at 0.6 AU from the central star. In the lower plot we show 
the intensity derived from our model on a log - log scale. The 
grey horizontal dashed lines marks the part of the red shifted 
part of the B2 spectrum that has been replaced as described in 
section |2] and the black horizontal dashed lines marks the part 
of the red shifted part of the B 1 spectrum of the removed tel- 
luric absorption line. 



drop in intensity as seen in HD101412. This disk has a flaring 
geometry, meaning that it has - depending on the disk models 
used - no or negligible shadowing by the inner rim. The flar- 
ing also allows the disk to intercept more stellar flux at larger 
distance, causing the [OI] emission up to a larger distance 
from the central star. The onset of emission for the B2 data 
of HD179218 is located beyond the dust sublimation radius. 
This result is rea ched by assuming the inclination found by 
Dent et al.1 d2005l) of 40+ 1 0° . When the lower limit of the inch- 
nation of 30° is assumed, the onset of emission starts at the dust 
sublimation radius. Both the shape and the inner plus outer radii 
of the [OI] emission are in reasonable agreement with those ex- 
pected from a flared disk model. The discrepancy between the 
observed and modeled emission at larg e radii is also seen in the 
two g roup I sources investigated by Ack e & van den Anckej 
(2006). This may be due to the model assumption that the OH 



abundance is homogeneous throughout the disk. 



as can be seen e.g. in the log - log part of Figure[T3] The the emission 
originating from 2 AU of the central star already has dropped in inten- 
sity by a factor of 100 compared to the maximum. 



Fig. 14. Same as Figure[l3]for HD179218 . The spikes at small 
radii and at 4 AU are an artifact of our data reduction and not 
significant. 

4.2.2. HD1 35344 B 

The spectrum of HD 135344 B is contaminated by very broad 
atmospheric absorption lines. Because of this and because of 
the fact that the [OI] feature is overlapped by two of these 
lines as shown in Figure Q3J we have undertaken one extra 
step to reduce these data and therefore will discuss this source 
separately. We correct the absorption by these two lines with 
two models: (1) A double Voigt model, because a Voigt pro- 
file closely mimics the expected profile of a stellar atmospheric 
line and (2) a polynomial model (a good fit was found using a 
6th order Chebyshev-polynomial), both shown in Figure Q2] 
A problem when fitting the Voigt profile was that we were 
unable to fit the absorption using the location of the lines 
found in HIP53070, a reference star with the same spectral 
type as HD135344 B (Figure [TBIl. Using these values results 
in over(under) fitting the right(left) wing of the absorption. For 
the further modeling of HD 135344 B we use these two mod- 
els and their average, as shown in Figure [17] This Figure also 
shows that the blue- and red shifted sides are asymmetric. This 
can either be caused by an asymmetry in the emitting region, 
or by an imprecise correction for the photospheric absorption 
(close inspection of Figure Q3] suggests that the blue side of the 
line lies on a monotonous rising slope, while the red side of 
the line may overlap with the minimum of the mixed absorp- 
tion line). Because of this, the blue- and red shifted sides are 
considered separately in Figures [18] and [19] There is a big dif- 
ference between the red and blue wing in the first AU, as well 
as between the different models and data procured at different 
times. The blue shifted profile shows a steep initial rise fol- 
lowed by a big drop resembling HD 101412, whereas the red 
side rises more slowly, and displays a smaller drop in inten- 
sity at a larger radius. There is also a big difference between 
the two models (and thus also their average) used in the first 
couple of AU. For the Chebyshev model the emission consis- 
tently starts at larger radii than that of the Voigt model. In all 
cases however, the emission extends up to a larger distance and 
with more brightness than the other two sources and vanishes 
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Fig. 15. The spectra of HD 135344 B (bottom) and HIP53070 
(top), another F4V star but with low v sin i. This plot shows that 
there are two broadened absorption lines marked with vertical 
lines mixing with the [OI] emission line of HD 135344 B . All 
the other absorption lines are caused by telluric absorption. 
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Fig. 16. The Chebyshev (dashed line) and the double Voigt 
(dot-dashed line) model used to model the photospheric ab- 
sorption in HD135344 B . The telluric absorption lines at the 
location of the four horizontal lines have been clipped out. 



around 90 AU. This value is comparable with t he outer radius 
of CO gas of 75 + 5 AU from bent et all d2005l) . 

Considering (1) that it is not possible to properly correct 
for the photospheric absorption, (2) the fact that HD 135344 B 
has an SED that is uncharacteristically for either a flaring or 
a self-shadowing disk because of the high (near) infrared ex- 
cess between 1 and 4 /mi and (3) the fact the peak of the [OI] 
emission velocity distribution does not coincide with the stellar 
photosphere we consider it likely that the observed [OI] emis- 
sion could also be influenced by other components than solely 
a disk such as an outflow with thermally exited oxygen. We 
will therefore limit our interpretation of our [OI] data on this 
star to a purely phenomenological description of the spectra. 



Fig. 17. The velocity corrected, normalized and converted to 
velocity around the [OI] line spectra corrected for the photo- 
spheric absorption using the 3 models of HD 135344 B . From 
top to bottom: the double Voigt model, the average model and 
the Chebyshev model. The horizontal bars mark the location of 
the clipped out telluric absorption lines. Only the spectra of Bl 
are shown 
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Fig. 18. The intensity versus radius graph for the blue (upper 
graph, triangles)- and red(lower graph, diamonds) shifted sides 
of HD 135 344 B showing the distribution of [OI] emission as a 
function of distance from the central star. The three lines rep- 
resent the Bl measurements for the Voigt model (upper grey 
line), the Chebyshev model (lower grey line) and the average 
of the two (thick black line), the light grey smooth line rep- 
resents a simple model for the [OI] emission from a disk in 
Keplerian rotation. 

5. Discussion and Conclusions 

We have detected non-thermal [OI] emission in all three tar- 
gets, one of which shows a double peaked emission profile 
and another one showing a a tentative double peak but cer- 
tainly flattened emission profile indicating emission from a cir- 
cumstellar disk. We see in the [OI] emission of one of these 
sources, HD101412, evidence for the existence of a puffed up 
inner rim followed by a shadow, supporting the disk model 
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Fig. 19. Same as Figure[l8]for the B2 measurements 



of iDullemond et al.l (12001b . Although this model can explain 
most of our observations there is one unexpected feature; the 
re-brightening (the 2 nd bump) of the disk of HD 10 14 12 shortly 
after the shadow cast by the inner rim. We expect the disk to re- 
main in the shadow much longer and possibly never emerging 
from it, thus preventing the disk to become more bright. We of- 
fer a possible explanation for this behavior by speculating that 
the dust could have a lower scale height than the gas, allow- 
ing for a layer of gas above the dust disk. This layer may well 
(partially) rise above the shadow of the inner rim causing the 
behavior seen in Figure [13] The [OI] emission of HD 1792 18 
is in good agreement with emission from the surface of a flar- 
ing circumstellar disk inclined at a 30 degree angle. Finally, we 
encounter several difficulties in explaining the [OI] emission 
of HD135344 B by a disk and therefore suggest the possibility 
of more components responsible for the emission than just the 
disk. 

We will improve on the modeling done in this paper in the 
next papers, where we will present MIDI observations of these 
stars showing warm dust enabling us to compare the coupling 
between gas and dust. We will then also model all our observa- 
tions together with the SEDs to further test and improve on the 
method employed in this paper, as well as to constrain as many 
stellar parameters on our targets as possible. 
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